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Abstract 



The assumption of exact, unbroken parity symmetry leads directly to a simple pre- 
dictive resolution of the atmospheric and solar neutrino puzzles. This is because the 
existence of this symmetry implies the existence of a set of mirror neutrinos which 
■ must mix maximally with the known neutrinos if neutrinos have mass. The maximal 

mixing of the electron neutrino with the mirror electron neutrino with 3 x 10~ 10 eV 2 ~ 
\5m 2 \ ~ 1CT 3 eV 2 leads to a predicted reduction of the solar neutrino flux by a factor of 
2, which is in quite good agreement with the experiments. The maximal mixing of the 
muon neutrino with the mirror muon neutrino with \5m 2 \ ~ 1CT 2 eV 2 also solves the 
atmospheric neutrino puzzle. We show that there is a significant range of parameters 
where these solutions are not in conflict with standard Big Bang Nucleosynthesis when 
the creation of lepton asymmetry due to neutrino oscillations is taken into account. 



It has been known for a long time[|IJ but not widely appreciated that it is possible to 
build a phenomenologically consistent gauge model which has a parity symmetry which 
need not be broken at all. In order to achieve unbroken parity symmetry it is necessary 
to double the number of fermions and the gauge symmetry. However, while the number 
of particles is doubled the number of parameters is not significantly increased (only 
two additional parameters in the minimal model with massless neutrinos)^]. If the 
neutrinos in the exact parity symmetric model have mass, and if mass mixing between 
ordinary and mirror neutrinos exists, then the mass eigenstate neutrinos must also be 
parity eigenstates. Because parity transformations simply interchange ordinary and 
mirror neutrinos, the mass eigenstate fields will be maximal combinations of ordinary 
and mirror weak eigenstates^]. This result holds independently of the details of the 
origin of the neutrino masses. 

The mirror neutrinos are essentially sterile as far as ordinary interactions are con- 
cerned. However unlike totally sterile neutrinos (such as right-handed neutrino gauge 
singlets) mirror neutrinos interact amongst themselves and with mirror particles with 
interactions of the same form and strength as ordinary neutrinos interact with ordi- 
nary particles. One advantage of mirror neutrinos over conventional sterile neutrinos 
is that the mirror gauge symmetry provides an excellent reason as to why they are not 
very heavy|2|, If the parity symmetry connecting the ordinary and mirror worlds is 
unbroken, then the mirror neutrinos are set by the same scale as ordinary neutrinos ||. 

We will denote the three mirror neutrinos by 

K, *v K- (!) 

With small intergenerational mixing it follows from the unbroken parity symmetry 
of the model that v e and v' e will be approximately maximal mixtures of two mass 
eigenstates v\ and v^- Similarly, z/ M and v' will each be approximately maximal mixtures 
of two mass eigenstates as will u T and v' T . We will denote the 5m 2 describing these 
maximal oscillations by 5m^,, Sm 2 , and 5m 2 T , respectively. 

The maximal mixing of v e and z/ will solve the solar neutrino problem for the large 
range of parameters f|, [5], || 

3 x 1(T 10 eV 2 £ \6m 2 ee ,\ £ 1(T 3 eV 2 . (2) 

The maximal mixing of the electron neutrino with the mirror electron neutrino with 
8m 2 e i in the above range leads to a predicted reduction of the solar neutrino flux by 
a factor of 2 which is in quite good agreement with the experiments. (See Ref. Q 
for a detailed comparison of the predictions of the exact parity model with the solar 
neutrino dataj7[). 

The deficit of atmospheric muon neutrinos can be explained if there are — v' 

oscillations with sin 2 29 ~ 0.5 and 10 -3 eV 2 ~ \5m 2 ^,\ ~ KT 1 eV 2 @, |. The best 
fit || occurs for sin 2 29o ~ 1 and 

\5m 2 \^lO~ 2 eV 2 . (3) 
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The exact parity symmetric model is also compatible with the LSND signal |TD|, 

A potential problem with any model that has additional light degrees of freedom 
is that these extra states can contribute to the energy density of the early Universe 
and spoil the reasonably successful Big Bang Nucleosynthesis (BBN) predictions. This 
presents a problem for the exact parity symmetric model because it can potentially 
lead to a doubling of the energy density at the time of nucleosynthesis (which is 
equivalent to about 6 additional neutrinos). However, it is plausible that an initial 
macroscopic asymmetry between ordinary and mirror matter might exist, as can be 
arranged through the inflationary scenario proposed in Ref . ||1 1|| (for example). Even if 
ordinary matter dominates mirror matter immediately after the Big Bang, the oscilla- 
tions between the ordinary and mirror neutrinos might be expected to bring the mirror 
sector into equilibrium with the ordinary particles For maximally mixed ordinary 
- mirror neutrinos, the following BBN bounds have been obtained |L3| assuming that 
the lepton number asymmetry could be neglected: 

\6m 2 ee ,\ £ 10" 8 eV 2 , \Sm 2 ^,\, \6m 2 TT ,\ £ 1(T 6 eV 2 . (4) 

With the parameter choices Eq.(fJ, |^), there is a potential conflict with the naive BBN 
bounds Eq.(|4]). However, these bounds do not hold if there is an appreciable lepton 
asymmetry in the early Universe for temperatures between 1—30 MeV0]. Remarkably, 
it turns out that ordinary-sterile neutrino oscillations can by themselves create lepton 
number |15 , |I~6[ |T7f . Recently, we have shown [pj]] that the lepton number generated by 



ordinary - sterile neutrino oscillations can allow the bounds in Eq.(|]) to be evaded by 
many orders of magnitude. Indeed, the bounds can be relaxed sufficiently so that the 
ordinary - sterile neutrino oscillation solutions to the solar and atmospheric neutrino 
anomalies do not significantly modify BBN. 

The purpose of this paper is to study the special case where the sterile neutrinos 
are mirror neutrinos 0. The mirror neutrinos are essentially sterile when probed by 
ordinary matter, however they do have significant self interactions. There are two main 
effects of the self interactions in the early Universe. First, the effective potential govern- 
ing ordinary - mirror neutrino oscillations will gain a contribution from the interactions 
of the mirror neutrino with the background. Second, the mirror weak interactions can 
bring the mirror neutrinos into equilibrium with the other mirror particles. This effect 
is quite important because it will significantly modify the momentum distribution and 
the number density of mirror neutrinos compared with the case of sterile neutrinos. In 
Ref. ]17|| , we showed in detail how lepton number creation can evade the naive bounds 



of Eq.(|^) and we determined the required parameter space for strictly sterile neutrinos. 
In the present work we consider the case of mirror neutrinos. For this case it turns 
out that the effect of the mirror neutrino self interactions is to significantly enlarge the 
allowed region of parameter space compared to the case of strictly sterile neutrinos. 

There are many independent 6m 2 , sin 2 29q parameters. We will need to make 
some assumptions otherwise we cannot say anything definite. We will assume that 
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intergeneration mixing is small and that 

m Ve , m u > e < m Ufi ,m^ < m UT , m^. (5) 

The assumption of small mixing between the generations is quite natural in our opinion 
in view of the situation with quarks. This assumption is also supported by the LSND 



experiment ||10||, which claims to have measured a small mixing between the muon and 
electron anti-neutrinos. The above assumption for the mass ranges of the neutrinos is 
also quite natural in view of the mass hierarchy between generations for the quarks and 
leptons. It is also compatible with the three experimental neutrino anomalies (solar, 
atmospheric and LSND). [Of course it is not the only possibility, but represents our 
best guess given the existing information]. 

For ordinary - sterile or ordinary - mirror neutrino two state mass mixing, the 
weak-eigenstates [y w z/ s ) will be linear combinations of two mass eigenstates {y ai v^)\ 

v a = cos 9 v a + sin 6*0^, v s = — sin d§v a + cos O^Vb- (6) 

Note we will always define 8 in such a way so that cos2# > 0. We also adopt the 
convention that 5m 2 = m 2 — m 2 a . Hence with this convention 5m 2 is positive (negative) 
provided that > m a (rrib < m a ). 

Ordinary - sterile or ordinary - mirror neutrino oscillations can generate significant 
lepton number in the early Universe |15|, |16[ |I7| . The origin of this phenomenon can be 
traced to the fact that the effective potential induced from the coherent forward scat- 
tering of neutrinos with the background is generally unequal to the effective potential 
for anti- neutrinos if the background is CP asymmetric (|l8l, [L9| . This means that the 
matter mixing angles for neutrinos are generally unequal to the matter mixing angles 
for anti-neutrinos, and thus the oscillation rates for neutrino oscillations need not be 
the same as the oscillation rates for anti-neutrino oscillations. 

We begin by briefly reviewing the case of ordinary - sterile neutrino oscillations as 
developed in Ref.|]T^|. The evolution of lepton number L„ a = {n Va — n^ a )/n 7 (where 



the n's are the number densities and a = e, /i, r) in the early Universe generated by 
v a — v s oscillations can be approximately described by the following equation |T7|], 



dL Ua 7i 2 r s 2 T p Ua aP(c - W){dn+ a - dn+) 
dt ~ 4C(3)T 3 J [ X p + (c -bP + aP) 2 ][xP +\c - fr- aP) 2 ] + ' U 

where A is a small correction term, 

a ~, f s 2 TP Ua [xP + (aP) 2 + (b p - c) 2 ](dn- a - dn~) 

8C(3)W [xp + (c - bP + aP) 2 }[xP + {c - bP - aP) 2 } ' 1 ' 

In these equations, C(3) is the Riemann zeta function of 3 (C(3) — 1.202), dn^ a = 
dn Ua ± dn Ua , and similarly for dn^ s . In the region where the lepton number is much 
less than 1, or equivalently, the chemical potential can be neglected, 

a a 1 P 2d P fa\ 
dn v ~ dn p ~ — t=. 9 

2vr 2 1 + eP/ T 
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In Eq.(|7], ^|), c = cos2# , s = sin 26*0, and the quantities, b P ,a p ,x p ,T P , are all functions 
of momentum of the form: 

rr'P — n 2 I r L ( V \ 2 -np _ pi _P_ 
X - 6 + 4{A^F V<P> J ' 1 "« ~ 1 <P> ' 

fl - A* > fa UUJ 

where n 1 = 2((3)T 3 /tt 2 , Aq = 5m 2 /2p and the thermally averaged collision frequencies 
Y Ua are 

r„ Q ~ y a G 2 F T\ (11) 

with y e ~ 4.0, y p : r — 2.90. In Eq.(|lO|), G F is the Fermi constant, M w is the W-boson 
mass and A e ~ 55.0, A^ :T ~ 15.3. The function L/ a ' is given by 

=L Ua + L Ue + L^+L UT + V , (12) 

where 77 is a small asymmetry term which arises from the asymmetries of baryons and 
electrons and is expected to be about 10~ 10 [|n|. With the definitions Eq.flTUD observe 
that Eq. (|7]) implies that significant lepton number can only be generated provided that 
Sm 2 < and in this case only for oscillations with V < c. 

In Ref . [fT7 1 , we showed that the distribution of sterile states is governed approxi- 
mately by the equations 



dz (1-z) Tlj 2 dz (1-z) Vis 2 



dt 4 xp + (c - bP + aP) 2 ' dt 4 xp + (c - bP - aP) 2 ' 

where 



(13) 



dn v Jdp _ _ dn PB /dp 



dn Ua /dp' ' dn 9 Jdp 
Thus, the equation governing the evolution of L Ua has the approximate form 

dL Va _ 1 s 2 T p Va oP{c-lP) (1 — z + )p 2 dp 

~dT ~ 4((3)T 3 Jo [xp + (c-bP + aP) 2 } [xp + (c - &p - aP) 2 } (1 + cp/ t ) ^ 

where A is a small correction term, 



_ 1 s 2 TP Va [xP + (aP) 2 + ^ - c) 2 ] gVrip 

8C(3)T 3 io [a;P + (c-6P + aP) 2 ][xf + (c-6f-aP) 2 ] (l + eP/ T )' 1 J 

with 2: = (2 ± z)/2 and we have neglected a small term proportional to L Va . 

In the case of ordinary - mirror neutrino oscillations, there are two important modifi- 
cations. First, the effective potential of the mirror neutrinos is generally non- negligible. 
Second, the mirror weak interactions can bring the mirror neutrino into thermal equi- 
librium with other light mirror particles. We now discuss these points in more detail. 

In the case of u a — 1A oscillations, the dynamics depend on the difference of the 
effective potentials, 

V = V a - V* p , (17) 
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where V a (Va) is the effective potential experienced by a pure weak (mirror) eigenstate. 
These effective potentials can be expressed in terms of the parameters a p , \f (a' p , b' p ) 
as follows, 

V a = (-a p + b p )A p , V' p = (-a' p + b' p )K (18) 

If the number of mirror neutrinos is much less than the number of ordinary neutrinos 
then b' p ~ 0. [Note that the 6-part of the effective potential is proportional to the 
number densities of the background particles. This dependence is not given explicitly in 
Eq. ( |I0D since for this equation the number densities were set equal to their equilibrium 
values]. The parameter a' p has the form 

(19) 

where L'^ is given by 

L'W = L v i + L v , + + L K + rf, (20) 

and L u > are the mirror lepton numbers, which are defined by L v i = {n v > — njj')/n 7 
(note that n 1 is the number density of ordinary photons) and rj' is a function of the 
mirror baryon/electron number asymmetries [which is defined analogously to 77]. We 
will assume that rf is small and can be approximately neglected. Thus, assuming that 
the number density of mirror particles is much less than the number density of ordinary 
particles, the modification of the effective potential due to the mirror interactions of 
v'q can be approximately taken into account by simply replacing in the definition 

of a p by L = L^ a > — L'^> . Since ordinary + mirror lepton number is conserved (and 
we will assume that it is zero), it follows that 

L Ve + L Vvt + L Vt + L K + L„, + L K = 0. (21) 

For example, if we consider u a — Up oscillations in isolation then = —L Va and 

L ~2L( Q ). 

Another important effect of the mirror interactions is that the momentum distribu- 
tion of the mirror neutrinos will be approximately Fermi-Dirac distributions, and the 
other mirror particles will be excited until 

T v ' e = T v i^ = T u ' r = T e i = Ty = T , (22) 

we consider the V ~ 100 MeV region only where we can approximately neglect the 
excitation of mirror muons and other heavier mirror particles). Of course the light 
mirror particles will only be excited provided that there are sufficient mirror neutrinos 
around so that the interaction rates will be faster than the expansion rate, that is, 

hG 2 F T' 5 ~ 5.5T 2 /M P , (23) 
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where Mp is the Planck mass and h is a numerical parameter which depends on the 
particular interaction (see Table 1 of Ref.|2(J f° r a ns ^ of the reaction rates). Typically, 
h ~ 1/8. Solving Eq. (B3T) , we obtain the condition 



T' ~ 2 



T 



MeV 



MeV. 



(24) 



Assuming that the above condition is approximately satisfied, the system of ordinary 
and mirror particles form two weakly coupled thermodynamic systems; the system 
comprising the ordinary particles at a temperature T, and the system comprising the 
mirror particles which has a distinct temperature T". 

Let us determine the equation governing the evolution of T'. Initially, we will 
assume that T' = 0. Ordinary - mirror neutrino oscillations can then generate a mirror 
neutrino v'q say. The rate at which mirror neutrinos are created / destroyed by u a — v'a 
oscillations is governed by the rate equation, 



d 
dt 



n v ' + n D > fj 



1 



n Ua + Hp, 

The reaction rates are given by]T7|, |22 

1 



T(v a -> v' p )(dn Ua - dn u + F(u a -> v' p )(dn Da - dn D 0. 

(25) 



r(v n 



-T* sin 2 2^ m (sin 2 ^— ) = - 

2 u a m\ 2L rnJ ± 



xp + (bP - aP - c) 5 



(26) 



The rate for for the process T(u Q 
the above equation. If Eq. 
neutrinos has the form 



dn 



v, 



can be obtained by replacing aP — > — a p in 
is satisfied, then the momentum distribution of sterile 

, ~ dn I P" d P (27) 

where we have neglected the mirror neutrino chemical potential which is small provided 
that the lepton number is sufficiently small. Thus using Eqs. (p6|) and Eq . fl27[) , it is 
straightforward to show that Eq.(B3|) can be expressed as follows, 



d_ 

dt 

where 



'n v ' + n p r 




f°° s 2 T p Ua [(b p - c) 2 + (a p f + x p ]p 2 F{p, T, T')dp 


n v + rip 


6C(3)T 3 J 


o [xp + (c-bP + aP) 2 } [x p + (c-bP - clp) 2 ] 



F(p,T,T') 



1 



(2f 



(29) 



1 + eP/ T 1 + eP/ T ' ' 

During the process whereby the mirror interactions excite the other mirror particles 
and thermalize the momentum distributions, the number density of sterile neutrinos 
is generally not conserved. Thus, we cannot directly use Eq.(^) to determine the 
evolution of X". However, the mirror interactions must conserve energy. If the energy 
density of the mirror particles is much less than the energy density of the ordinary 
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particles then the process whereby the mirror sector comes into equilibrium with itself 
should not significantly affect the expansion rate of the Universe. For this reason, 
and because of the conservation of energy, it follows that the energy density of the 
mirror particles normalized to the energy density of the ordinary particles should to 
a good approximation not change due to the expansion of the Universe. We will 
denote this quantity by 7 = p'/p. In the region where 1 MeV ~ T ~ 100 MeV, and 
1 MeV ~ T' ~ 100 MeV, then 7 ~ j"4^y4 (assuming that we are in a region where 



Eq. (p4|) is valid). Using Eq . (|28|) , it is straightforward to show that the evolution of 7, 
which is entirely due to ordinary - mirror transitions, satisfies the following equation 
(where we are assuming that 7 <C 1), 

dl _ 1 f°° ^ K 6 " - c ) 2 + ( gP ) 2 + x P ]P 3 HP, T, T')dp 



dt 19C(3)A^T 4 7o [xp + (c - 1? + aP) 2 ][xP + (c - & - aP) 2 ] ' 

where N = p' / p' v ,. Using p' = 3p' u , + p' e , + p' , = §\p' v ,, it follows that N ~ 6.14. 

In the region where Eq . (^4[) holds, Eq.(Q) and Eq.([27|) imply that the evolution of 
lepton number obeys the following equation, 

dL Va 1 [°° s 2 TP /a gP(c-bP)p 2 F(p,T,r)dp 

dt ~4C(3)T 3 7o [xp + (c - bP + aP) 2 }[ X P + {c - bP - aP) 2 ] ' 1 ' 



where [23 



and 



2tt 2 Jo [xp + (c-bP + aP) 2 ] [xp + (c - \P - aP) 2 \ ' 1 ' 



1 p 2 e p/T 1 p 2 e p/r 

G(p, T, T ) = — (i - ep/T)2 - — p - ep/r/)2 . (33) 



Equations (|30| ) and are coupled differential equations that can be solved for L Ua , V . 



T'<^2(-f—) MeV, (34) 



In the region where 

T' <f 9 , 

VMeV, 

the number of mirror particles are insufficient to enable the mirror sector to come into 
equilibrium with itself. In the case of Eq. (|34|) , the momentum distribution and number 
density of mirror neutrinos are not significantly affected by the mirror interactions. This 
means that the momentum distribution of mirror neutrinos should be the same as with 
the case of sterile neutrinos. Thus the momentum distribution and number density 
of mirror neutrinos should be governed approximately by Eq.([T3|) with the evolution 
of lepton number governed approximately by Eq. ([T5|) . Of course even in the case of 
Eq . (|34"D , the effects of the mirror interactions on the effective potential must still be 

taken into account, that is, we must replace L (a) by L 
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There are several ways in which the creation of lepton number (s) can prevent ster- 
ile/mirror neutrinos from coming into equilibrium. One way is that one set of oscilla- 
tions u a — v s creates L Ua . The lepton number L Va can then suppress other, independent 
oscillations such as Vp — v s oscillations (with /3 ^ a) for example. A more direct, but 
less dramatic way in which the creation of lepton number can help prevent the ster- 
ile/mirror neutrinos from coming into equilibrium, is that the lepton number generated 
from say u a — v s oscillations itself suppresses the v a — v s oscillations |Q. We will ex- 
amine the latter effect here (the former effect will be studied in a moment). Previous 
work[|l3j obtained the BBN bound on ordinary - sterile neutrino oscillations with large 
\5m 2 \ ~ 10~ 4 eV 2 (with 5m 2 < 0) and small sin 2 26q ~ 1CT 2 . This bound can be 
approximately parametrized as follows 

(35) 

This bound arises by assuming that the v a — v s oscillations do not bring the sterile v B 
state into equilibrium. Note that this bound neglected the creation of lepton number 
and it also did not include the effects of the distribution of neutrino momenta. However, 
in the realistic case, the creation of L Va (after it occurs), will suppress the v a — v s 
oscillations and the actual bound might be expected to be somewhat less stringent 




than Eq. (p5|) . Nevertheless, for the case of truly sterile neutrinos we found |L7| that 
Eq. fl3"5"D turned out to be a reasonable approximation in the realistic case where the 
momentum distribution and the creation of lepton number was taken into account. This 
is largely due to the fact that the creation of sterile neutrinos suppresses the lepton 
number creating oscillations and consequently delays the point where significant lepton 



number can be created [Oj. In the case of mirror neutrinos, there will be much less 



of this type of effect because the number density of the mirror neutrinos is kept low 
due to the excitation of the other mirror particles. Also, the thermalization of the 
neutrino momentum distributions means that not all of the mirror neutrinos will have 
low momentum. Thus, we would expect that the bound Eq . (|35|) should be weakened 
somewhat. 

For the case of ordinary - mirror neutrino oscillations, the BBN bound for large 
\Sm 2 \ ~ 1(T 4 eV 2 (with Sm 2 < 0) and small sin 2 28 ~ 1(T 2 can be obtained by solving 
the coupled differential equations, Eq . (|3~0D and Eq . fl3T|) . Doing this numerical exercise, 
we obtain the following bounds assuming that 8N e ff = p' / p Ua < 0.6 (1.5) |2T]] (where 
N e ff is the effective number of neutrinos present during nucleosynthesis), 

sin 2 26*o ~ 3(6) x 10 -4 (]f^j) 2 f° r v h,t ~~ v 'p oscillations, 
sin 2 26 ~ 1(2) x 10~ 4 (jg^) 5 for v e - v' p oscillations. (36) 



Comparing these bounds with Eq. (B5p , we see that the BBN bound on sin 2Qq for 
ordinary - mirror neutrino oscillations is about an order of magnitude weaker than the 
corresponding bound for ordinary - sterile neutrino oscillations. 
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We now identify the region of parameter space for which the maximal ordinary - 
mirror neutrino oscillations can solve the solar and atmospheric neutrino anomalies 
without leading to a significant modification of BBN. 

We first study the maximal ordinary - mirror neutrino oscillation solution to the 
atmospheric neutrino problem. We will assume that the various oscillations can be 
approximately broken up into the pairwise oscillations u a — v*. We will denote the 
various oscillation parameters in a self-evident notation, 

^aB'i a a/3' f° r v a — v'a oscillations, (37) 

where a, (3 — e,fi,r. We will denote the mixing parameters, 6m 2 , sin 2 2#o appropri- 
ate for v a — Up oscillations by bm 2 a ^,, sin 2 20Q 13 . Note that lepton number cannot be 

created by u a — v'^ oscillations until (kP a ai) ~ cos26'q^ . Recall that the V parameter 
is inversely proportional to 5m 2 [see Eq.(|II])]. Thus, the earliest point during the evo- 
lution of the Universe where lepton number can be created due to ordinary - mirror 
neutrino oscillations occurs for oscillations which have the largest |<5m 2 |. Note that 
these oscillations should satisfy the BBN bound Eq.(PB[) if we demand that the sterile 
neutrino should not significantly modify BBN. Note that the — v' oscillations have 

quite small |5m 2 M ,| ~ 1(T 2 eV 2 , and cos 29^ ~ (assuming; maximal or near maxi- 
mal mixing), and thus these oscillations themselves cannot produce significant lepton 
number. However, the \Sm 2 \ for v T — v' (or v T — z/) oscillations can be much larger. 
Also note that 5m 2 < if m Vr > m v i (or m Vr > m u ij. 

We will first consider the system comprising z/ T , z/ M and u' (and their anti-particles). 
The effects of the other neutrinos will be discussed in a moment. We will assume for 
definiteness that m„ T > m v i so that l^m 2 ^! > l^m 2 ^,] and the v T — i/ oscillations 
create L Vt first (with L Ufi ,L Ue assumed to be initially negligible). 

Recall that in the case of v a — v'^ oscillations, the effect of the mirror interactions 
on the effective potential can be taken into account by replacing in the definition 

of oP by L = L^ a > — L'^'. For the v T , z/ M , v' system, 

L ~ 2L, M + L Vt - 2L^ ~ 4L, M + 3L, T , 

L Tfl ^ 2L Ut + L Ufi - 2L v , t ~ AL Ut + 3L, M , (38) 

where we have used Eq.(^T|) with L Ve ~ L^/ ~ — 0. Thus, the creation of L Ut 

by u T — u' oscillations also implies the creation of the quantity L ■ Because aP^, is 

directly proportional to L it follows that the creation of a large L will make a? , 
large (i.e. significantly greater than 1), which thereby suppresses the oscillations (note 
that sin 2 29^ ~ sin 2 29^ /a 2 ^/ if a^i ^> 1). However, note that — u 1 oscillations 

can potentially destroy L (because for these oscillations L — is an approximately 
stable fixed point for the temperature range of interest). Thus, we must obtain the 
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j nit 



region of parameter space where L created by v T — v' oscillations does not get sub- 
sequently destroyed by — v' oscillations. We will determine this region of parameter 
space by numerically solving the coupled differential equations governing the evolution 
of L Ur ,L Ufi and T' . [An approximate analytic computation could be done following 
similar reasoning to the case of sterile neutrinos ||17||. Note however that the approxi- 



mations used in this analytic derivation are not valid if the point where the destruction 
of L due to — 1/ oscillations reaches a maximum can occur during the phase where 

the creation of L due to v T — v' oscillations is still growing exponentially]. 

The rate of change of and L Ut due to the v T — u f , — u' oscillations can be 
obtained from Eq.pl]). This leads to the following coupled differential equations, 



dL^ j_ fOQ S ^2e^Tl^{cos2ef~bl^F(p,T,T')dp 

dt — 4C(3)T3 JO , x p + ( cos2 8^'-b p ,+a p ,) 2 ][x p ,+(cos 26?/ ~b p ,-a p ,) 2 ] 

L MM v U l^fJ' MM u MM MM 

_3L^ f00 sin 2 2^'r^[^ ; +(a^) 2 +(^,-cos2^') 2 ]G(p ) T,T')* 
2t 2 J ° [x p . + (00829^'-^ ,+a p .)^][x p . + (cos 29^' - b p ,-a p ,) 2 ]' 

L MM u MM MM MM u MM MM 

dL Ur 1_ f oo sin 2 20f V P T a p ^, (cos 29 f -b p ^, )p*F{p,T,T>)dp 

dt — 4C(3)T^ JO r P + i cos2 9 Tp ' -b p ,+a p ,) 2 }[x p ,+(cos26>^'-6 p ,-a p ,) 2 1 

tm u tm tm tm u tm tm 

_3L^ f oo sin 2 2 e -'r^[< fi/+ (a^ ; ) 2 + (6^-cos 2e - , ) 2 ]G( P ,T,T')d P 

2 *- 2 J ° [z p , + (cos2^'-6 p ,+a p ,)2][x p , + (008 20^'-^ ,-a p ,) 2 ]' ^ ' 

Tfl 1 V (J TfJ, TH f/ J L TM U TM TM 

where the equation governing the evolution of T" can be obtained from Eq.(^), and it 



is 



^ 1 ,oc sin 2 2^'^ [(Q -cos 2C') 2 + (<„, ) 2 +< ( ,,] P 3 F(p,T,r')<i P 

19C(3)7VT 4 JO r p + ( cos 26» f ^'-f,f ,+a p ,) 2 ][x p , + (cos 2^' - 6 P ,-a p ,) 2 ] 

L MM u MM MM MM MM MM 

i f00 sm 2 2 ^ , rg T [(^ M ^cos2e ( 7 , ) 2 + « M/ ) 2 +< M/ ] p 3 F(PiT;r)dp 

19C(3)AT 4 JO , x p + / cos29 ;/-b p ,+a p ,)' 2 }[x p ,+(cos29j/ -b p ,-a p ,) 2 ] ' ^ ' 

tm u tm tm tm u tm tm 



Note that Eq.(fMj) and Eq . fl39|) are coupled differential equations which must be solved 
simultaneously. In our numerical integration of Eqs. (|39|, [40]) , we will assume for def- 



initeness that sin 2 29^ and that 5m/ 1 , are given by the best fit for the atmospheric 
neutrino data, i.e. sin 2 26^ ~ 1 and {Sm 2 ^,] ~ 10~ 2 eV 2 . The results of this exercise 



are shown in Figure 1 |25| , pq] . In the region above the solid line, the L created by 
v T — i/ oscillations does not get subsequently destroyed by — i/ oscillations. The 
numerical work also demonstrates that the lepton number is generated early enough 
and is large enough to suppress the — v' oscillations sufficiently so that these oscil- 
lations do not create a significant number of i/ states. The requirement that v T — v' 

oscillations do not produce too many sterile states implies an upper limit on sin 2 2#q M 

[see Eq.(^)]. This upper limit has been shown in the Figure assuming for definiteness 

that p' / p Ua < 0.6 (dashed-dotted line). Also shown in Figure 1 (dashed line) is the 

,2 i < 1 «nn ^T/2n 



cosmological energy density bound \5m 2 Tll , \ ~ 1600 eV 2 §27§. 
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Comparing the allowed region of parameter space shown in Figure 1 with the analo- 



gous case for sterile neutrinosJT7j, it is apparent that the corresponding allowed region 
for mirror neutrinos is somewhat larger than the allowed region for sterile neutrinos. 
The increase of parameter space (which is about an order of magnitude larger for 
sin 2 2#o M ) in the case of mirror neutrinos is primarily due to the result that the bound, 
Eq.(|36D, is considerably less stringent than the bound, Eq . (p5|) . 



While we have focussed on the L generated by v T — v' oscillations, similar results 
will also hold for L generated by v T — v e oscillations. Note that in this case, 

L ~2L^+L Ut -2L^-L k ~4L^+2L„ t , 

L e ~ 2L Ut + L„„ - 2L„, - ~ AL Vt + 2L^. (41) 

In this case, solving the appropriate coupled differential equations, we find that the 

L generated by v T — v' e oscillations does not get destroyed by — u' oscillations for 
the range of parameter space shown in Figure 2. Of course we only require that the 
oscillation parameters are in the allowed region shown in Figure 1 or Figure 2. 

Note that our results should only be considered as approximate because we have 
only included the v T — z/g and — v' oscillations (where (3 = \x for Figure 1 and (3 = e 
for Figure 2). In the most general case, the system is considerably more complicated. 
In general one should also include the effects of the oscillations with 5m? > as well 
as ordinary - ordinary and mirror - mirror neutrino oscillations. The oscillations with 
5m? > 0, include oscillations such as — v' T [given the assumption Eq.(^)]. Focusing 
on this example, these oscillations generate L v ,L v i such that — L'^ — ► 0. They 

cannot prevent v T — v or v T — v e oscillations from generating L Ut and hence L ■ For 
this reason these oscillations cannot change anything qualitatively and for simplicity 
we have neglected them. The effect of ordinary - ordinary neutrino oscillations, such as 
u T — oscillations is to generate L Vji , L Vt such that L v — L Vt — > 0. These oscillations 

also cannot prevent L being generated by v T — v (or v T — v e ) oscillations. Also the 
rate of change lepton number due to these oscillations is typically suppressed compared 
to the rate of change of lepton number due to ordinary - sterile (or mirror) neutrino 
oscillations |fr5|, |T7| |. Finally, the effect of mirror - mirror neutrino oscillations may also 



be important, although the precise effect of these oscillations is less clear. 

We turn to a brief discussion of the maximal oscillation solution to the solar neutrino 
problem. Recall that maximal v e — v' e oscillations can lead to a simple predictive solution 
to the solar neutrino problem for the large range of parameters, Eq.(^). Note that much 
of this parameter space is naively in conflict with BBN [see Eq.([|)]. However, v T — u' e 

or v T — u' oscillations will generate L which can suppress the v e — u' e oscillations 

~ee' 

provided that u e — u' e oscillations do not destroy the L asymmetry. The situation is 
completely analogous to the case involving the — u' oscillations that we have been 
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studying above. For such large |5m^„/| or |5m^ e ,| identified in the Figures, it turns 
out that essentially the entire parameter space Eq.(@) does not lead to any significant 
modification of BBN (which is largely due to the fact that |5rn,g e /| ~ \Sm 2 ,\). 

In summary, the solutions of the atmospheric and solar neutrino problems suggested 
by the exact parity model are not in conflict with BBN for a significant range of pa- 
rameters. Consistency with BBN does require relatively large values of the parameters, 

\5m 2 ,\ (or \5m 2 Te ,\) ~ 25 eV 2 (50 eV 2 ). (42) 



This suggests that m Vr ~ 5 eV. Such values of the tau mass can also be motivated 
independently by the evidence for dark matter. Thus, demanding that the exact parity 
model solution of the atmospheric neutrino anomaly be consistent with BBN implies 
that the tau neutrino mass should be in a cosmologically interesting range. Also 
note that this tau neutrino mass range will be probed by the NOMAD and CHORUS 



experiments 28 . 



Figure Captions 

Figure 1. Region of parameter space in the sin 26^ , —drn 2 .^, plane where the L 
created by v T — v' oscillations does not get destroyed by — v' oscillations. This region 
which in the Figure is denoted by the 'Allowed Region" , includes all of the parameter 
space above the solid line. We have assumed that sin 2 29^ ~ 1 and l^w^ /| = 10~ 2 eV 2 
(which is the best fit to the atmospheric neutrino data). Also shown (dashed line) is the 
cosmology bound {Sin 2 .^, | ~ 1600 eV 2 . The dashed-dotted line is the nucleosynthesis 
bound Eq.(p6|) assuming that N e ff = p' '/ ' p Ua ~ 0.6. 

Figure 2. Region of parameter space in the sin 26^ , —8m 2 e , plane where the L 
created by v T — z/ oscillations does not get destroyed by — v 1 oscillations. As in 

Figure 1, we have assumed that sin 2 26*0^ ~ 1 and |5?n 2 M /| = 10~ 2 eV 2 . Also shown 

(dashed line) is the cosmology bound {Sm 2 .^ \ ~ 1600 eV 2 . The dashed-dotted line is 

the nucleosynthesis bound Eq.(|36"D assuming that 5N e ff = p' / p Ua ~ 0.6. 
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